Negative refractive index response of weakly and strongly coupled optical 

metamaterials. 
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We present a detailed study of the retrieved optical parameters, electrical permittivity, e, magnetic 
permeability, /x, and refractive index, n, of the coupled fishnet metamaterial structures as a function 
of the separation between layers. For the weak coupling case, the retrieved parameters are very close 
to the one-functional-layer results and converge relatively fast. For the strong coupling case, the 
retrieved parameters are completely different than the one unit fishnet results. We also demonstrate 
that the high value of the figure of merit (FOM = |Re(n)/Im(n)|) for the strongly coupled structures 
is due to the fact that the real part of the negative n moves away from the maximum of the imaginary 
part of n (close to the resonance), where the losses are high. 



I. INTRODUCTION 



Metamaterials are artificially engineered structures 
that have properties, such as negative refractive index, 
n, nonexistent in natural materials. The recent develop- 
ment of metamaterials 0] with negative n confirms that 
structures can be fabricated and interpreted as having 
both a negative permittivity, e, and a negative perme- 
ability, /i, simultaneously. Since the original microwave 
experiments for the demonstration of negative index be- 
havior in split ring resonators (SRRs) and wire struc- 
tures, new designs have been introduced, such as fishnet, 
that have pushed the existence of the negative refraction 
at Thz and optical wavelengths 0,Bi,IS[S0- Most of 
the experiments with the fishnet structure measure trans- 
mission, T, and reflection, R, and use the retrieval proce- 
dure 0,11 GHHH] to obtain the effective parameters, e , /x , 
and n. Although, stacking of three [10], four [l3j], 10 [ll| 
functional layers and recently fabricated 10-functional 
layer fishnets (21 layers of silver and MgF 2 ) have been re- 
alized, they do not constitute a bulk metamaterial. Even 
the thickest fabricated fishnet structure only has a 
total thickness, 830 nm, half of the wavelength (A=1700 
nm). Here, we report a detailed study of the weakly 
and strongly coupled fishnets to understand the origin of 
negative n, as well the mechanism of low losses (that is, 
high figure of merit (FOM)) for the weakly and strongly 
coupled fishnets. We also study the convergence of the re- 
trieval parameter (e, /x and n) as the number of unit cells 
(layers) increases. For the weakly coupled structures, the 
convergence results for n and FOM are close to the single 
unit cell. As expected, for the strongly coupled struc- 
tures, hybridization is observed and the retrieval results 
for n and FOM are completely different from the single 
unit cell. We demonstrate that the high value of FOM 
for the strongly coupled structure is due to the fact that 
the real part of negative n moves away from the maxi- 
mum of the imaginary part of n (close to the resonance) , 
where the losses are high. 



The idea of left-handed materials, i.e., materials with 
both negative e and negative /x, where the electric field 
(E) , magnetic field (H) , and wave vector (k) form a left- 
handed coordinate system was developed by Veselago [lj| 
decades ago. However, it was only recently that such ma- 
terials were investigated experimentally at high frequen- 
cies 0, H, |3, H, IS 0], an d the field is driven by a wide 
range of new applications, such as ultrahigh-resolution 
imaging system , cloaking devices [l7], Ell , and quan- 
tum levitation [19| . Realizing these applications, several 
goals must be achieved: three-dimensional rather than 
planar structure, isotropic design, and reduction of loss. 

Most of the metamaterials exhibiting artificial mag- 
netism [H, [13, US HE Hl| and a negative refractive index, 
n, at THz and optical frequencies [S H, 0, H, HI], consist 
of only a functional layer. The number of actual layers 
M = 2 x N+ 1, where N is the number of functional lay- 
ers. The first five-functional-layer of SRRs operating at 6 
THz was published 11411 in 2005, and four layers of SRRs 
operating at 70 THz [l3| was published in 2008. The first 
three-functional-layer of fishnets (7 layers of silver and 
MgF 2 ) operating at 200 THz was published [Hj] in 2007, 
and recently a 10-functional-layer of fishnets (21 layers 
of silver and MgF 2 ) operating at 200 THz was fabricated 
[7J. However, it is very important to study how the op- 
tical properties (e, \i and n) change as one increases the 
number of layers. How many layers are needed to achieve 
convergence of the optical properties and one can call this 
metamaterial bulk? How do optical properties behave as 
one changes the distance between two neighboring fish- 
nets? If the distance is small, we have a strong coupling 
case, and one achieves the photonic crystal limits. The 
convergence of optical properties is slow, and more im- 
portantly, it does not convergence to the isolated fishnet 
case. What is the mechanism for negative n in the strong 
coupling limit? 

In this paper, we present a detailed study of the re- 
trieved optical parameters, e, /x, and n of the single fish- 
net metamaterial structures as a function of the size of 
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the unit cell. We find that as the size of the unit cell de- 
creases, the magnitude of the retrieved effective param- 
eters increases. In order to understand the underlying 
physics of the coupled structures, we study the retrieved 
parameters of the coupled fishnets as a function of the 
distance between them. Finally, we study the conver- 
gence of the retrieved parameters as the number of the 
unit cell increases for the weakly and strongly coupled 
structures. For the weakly coupling case, the retrieved 
parameters are very close to the one-functional layer re- 
sults and converge relatively fast. For the strong coupling 
case, the retrieved parameters are completely different 
than the one unit fishnet results. The strong coupling 
case explains the recently observed negative refractive 
index in the 21-layer fishnet structure p}, especially the 
high FOM, due to the periodicity effects, as will be shown 
below. 



II. WEAKLY AND STRONGLY COUPLED 
FISHNETS 




FIG. 1: (a) Schematic of a fishnet structure with 11 metallic 
layers, (b) a single unit cell with geometric parameters marked 
on it. 

In Fig. [Tjwe present a schematic graph of the unit-cell 
of the fishnet structure. The size of the unit cell along 
the propagation direction is a z . a z is larger than the sum 
of the thickness of the metallic and the dielectric layers 
2t + s, where t and s are the thicknesses of the metal 
and dielectric layers, respectively. Notice the propaga- 
tion direction is perpendicular to the plane of the fishnet. 

In most of the experiments measuring the T and R of 
the fishnet structure [1, H, 0, H, d, H2| , there is only one 
layer of the sample measured. In this case, the unit cell 
size along the prop agation direction, a z , is undefined. 
We have shown [23|] that, as a z decreases, the magnitude 
of the retrieved parameters increases. It is well known 
from electronic systems that a monolayer of a surface 
can exhibit different properties from the bulk (many 
layers). So it is very important to systematically study 
whether the optical parameters of a single layer really 
correspond to the many layers system. We will study the 



weak and strong coupling limit of the two-layer fishnet 
structure. 



1 1 ' 1 1 1 ' 1 ■ 1 ' 1 1 










«\ f^iC f '•' ' ' 
\/\ V 

/ » \ / > i ■ 
\ / » 1 > 1 








a z 


\ / > ' ' ' 

\ / \ '. 
^ \ / \ ■ 

1 V/ • ' 










■ ! 

1 i \ i 




1 I I I I I I I 1 I I I I 

1.98 1.99 2.00 2.01 2.02 2.03 2.04 2.05 



FIG. 2: Retrieved real part of refractive index, n, from 
simulated data using unit cell size in the propagation direc- 
tion a z = a 1 15 (red), a z = 2a/ 15 (green) and a z = 4a/ 15 
(blue). Both one layer (dashed) and two layers (solid) re- 
sults are shown. The distances between two unit cells are 
d = a z — (2t + s) =0.04a, 0.11a, and 0.24a, respectively. 
The other geometric parameters are given by a x — a y — a, 
w x = 4a/15, w y — 3a/5, s = a/60, t = a/300, and the dielec- 
tric constant of the spacer is e r = 5. 

Figure [2] shows the real part of the effective refractive 
index, Re(n), as a function of A/a, for one layer and 
two layers of the fishnet structure described in Fig. [TJ 
with different distances between the unit cells. Notice 
the normalized resonance wavelength \ m /a 2.02, i.e. 
wavelength with maximum |Re(n)|, for one layer shifts 
only slightly when the size of the unit cell increases, but 
the magnitude of |Re(n)| decreases dramatically. For the 
two layers, when the distance, d, between them is large 
(d/a=0.2A, blue solid curve), the coupling between the 
two layers is weak and, therefore, the refractive index, 
Re(n), approaches the one layer simulation results. 
When the distance between the two layers becomes 
smaller (d/a=0.04, red solid curve) and the coupling- 
becomes stronger, hybridization takes place and two 
resonance modes exist, one at A/2 = 2.005, which gives 
Re(n) < 0; and one at 2.040, which has Re(n) > 0. 
The difference in value of the two resonance frequencies 
becomes larger as the distance between them decreases. 
Another very important issue is how fast the optical 
retrieval properties (e, /i and n) converge as the number 
of unit cells increases. We will present results for two 
cases, one for the weakly coupled fishnets. 

The only design that gave negative n at THz and 
optical frequencies is the so-called "double-fishnet" 
structure, which consists of a pair of metal fishnets 
separated by a dielectric spacer 0, H, 0, [E B 13 ■ F° r the 
incident polarization shown in Fig. [TJ the thin metallic 
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wires along the x-axis, parallel to the incident electric 
field, E, excite the plasmonic response and produce 
negative permittivity e up to the plasma frequency. 
Negative /i is obtained from the wires along the y-axis, 
parallel to the incident magnetic field H. At the 
magnetic resonance frequency, the two parallel bars 
sustain anti-parallel currents (along x-axis), providing 
a magnetic field B', mainly between the plates and 
directly opposite to the external magnetic field, H. The 
electric field, because of the opposite charges accumulate 
at the ends of the two metallic bars, is expected to be 
confined within the space between the plates and near 
the end points. Indeed, obtained simulations confirm 
this picture. 
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FIG. 3: Retrieved real part of effective refractive index, Re(n) 
for one layer (red solid), four layers (blue dashed), eight layers 
(green dotted) and ten layers (black dash-dotted) of the fish- 
net structure. The geometric parameters are a x — a y — 860 
nm, w x = 565 nm, w y — 265nm, s = 50 nm, t = 30 nm, 
d = 90 nm, and the spacer is made from MgF 2 with the di- 
electric constant e r = 1.9. The functional layers are separated 
by vacuum layers with thickness do as shown in the inset. 

In Fig. [3] we present the retrieved results for the ef- 
fective refractive index, Re(n), as a function of A for dif- 
ferent numbers of functional layers (N=l, 2, 3, 4 and 5) 
for weakly coupled fishnets system. The parameters are 
exactly the same as the strongly coupled case, that will 
be discussed below, but the spacing between the func- 
tional layers is d = 90 nm. As can be seen in Fig. [3l the 
retrieved results for Re(n) converge very fast (N=2) and 
the convergence results agree with the results of the one 
functional layer of the fishnet. 

When the fishnets strongly interact, it's not clear 
what the mechanism is for giving negative n. As 
discussed in Fig. [21 the isolated fishnet resonance 
frequency hybridizes into two different modes. The 
antisymmetric mode gives weak resonance with n w 0, 
while the symmetric mode gives a strong resonance 
with a strong negative n. In Fig. [U we present results 
for the retrieved, Re(n), for different number of layers 
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FIG. 4: The retrieved real part of n for 3, 7, 11, 19 and 
27 layers strongly coupled fishnet structure. The geometric 
parameters are a x = a y = 860 nm, w x — 565 nm, w y = 
265nm, s — 50 nm and t = 30 nm, and the spacer is made 
from MgF 2 with the dielectric constant e r — 1.9. The shadow 
region shows where the discontinuity happens. 



(3 to 27) for the recently fabricated Q negative index 
structure. Notice in the low wavelength limit (between 
1200 to 2100 nm), convergence of n is obtained and 
agrees with experimental results of Ref. 0] • In the high 
wavelength limit (A > 2200 nm), the Re(n) is zero and 
the Im(n) is much larger than the Re(n), exhibiting 
metallic behavior and transmission is equal to zero. This 
metallic behavior can be also seen in the transmission, 
T, (see the supplementary material) for the many layer 
structure. Above 2200nm, T is low, and behaves as a 
metal, while for A < 2000 nm, T is relatively large (~ 
0.8) and has Fabry-Perot resonances structure. The 
|Re(n) shown in Fig. Q] converges between 1200 nm and 
2200 nm to a finite value (positive for wavelengths less 
than 1500 nm and negative for 1500 nm < 2200 nm). 
For A >2200 nm, the |Re(n)| is zero and the |Im(n)| 
is large of the order of 3, and as expected for large 
wavelengths this strongly coupled metamaterial behaves 
as a metal. In addition, in Fig. [4] the 3-layer structure 
(the single fishnet structure) gives results completely 
different than those for the strongly coupled fishnets. 
These single fishnet results agree with those presented 
in Fig. [3] Another important quantity is the figure 
of merit (FOM) which can be defined two different 
ways. The usual definition is FOM=|Re(n)/Im(n)| 
and the experimental definition of Im(n) is given by 
Im(ra) = (A/4trZ) ln[(l - \R\)/\T\], where A, d, R, and T 
are the wavelength, sample thickness, reflectance, and 
transmittance, respectively. 
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FIG. 5: The figure of merit (FOM) of Re(n) < region for 
3, 7, 11, 19, and 27 layers strongly coupled fishnet structure. 
The FOM is calculated by FOM=|Re(n)/Im(n)|, where Re(n) 
is obtained by a retrieval procedure and Im(n) is calculated 
by Im(ra) = (A/47rd) ln[(l - \R\)/\T\}. 



III. FIGURE OF MERIT CALCULATIONS 

In Fig. \5\ we present the results of the FOM as a 
function of wavelength for different number of layers. 
For the one unit cell fishnet (3 layers), the FOM is really 
small (of the order of 2) and is located at A = 2100nm, 
the resonance frequency of the single fishnet structure. 
As the number of layers increases, the FOM increases 
and finally saturates to a constant value of the order of 
10. This behavior of the FOM for the strongly coupled 
fishnets is completely different for the weakly coupled 
fishnets, where the FOM does not change dramatically 
[23| as one uses more unit cells. Why is the FOM in 
the strongly coupled fishnets so much different than 
the single fishnet? It has been argued @, H3] that the 
FOM is larger because of the strong coupling between 
the neighboring layers, which provides destructive inter- 
ference of the antisymmetric currents across the metal 
film and effectively cancels the current in the center of 
the film, and, therefore, reduces the losses. We have 
systematically studied the current density for the differ- 
ent number of strongly coupled fishnet structures. For 
the single fishnet structure, the current density is along 
opposite directions in the two metallic bars. This is the 
typical behavior of negative index materials. When the 
number of layers increases, the current density is more 
complicated and there is no clear physical explanation 
why one obtains negative n and why the FOM is so large. 

In Fig. [5J we present the current density along the 
x-axis (or E-direction as shown in Fig. [I}, J x , of the 
antisymmetric and symmetric modes for the seven layers 
(4 metallic layers and 3 dielectric layers) strongly cou- 
pled structure. For the antisymmetric mode (as shown 
in Fig. [S](a)), two double- fishnets are formed by the first 



FIG. 6: (a) The current density distribution for a 7 lay- 
ers strongly coupled fishnet at wavelength, A = 2230 nm 
(antisymmetric mode), with Re(n) = —0.17. (b) The cur- 
rent density distribution for a 7 layers strongly coupled fish- 
net at wavelength, A = 1859 nm (symmetric mode), with 
Re(n) = —2.5. The cross-section is perpendicular to the y- 
axis (i.e., incident magnetic filed, H, direction). The color 
shows the current density in x-direction, J x , with the red and 
blue being the positive maximum and negative maximum of 
J x , respectively. The arrows show the direction of current 
density inside the silver layers schematically. 



and second silver layers, and by the third and fourth 
silver layers. The induced current inside two double- 
fishnets excite the magnetic fields, B', along the same 
direction. However, the second and the third silver lay- 
ers also form a double-fishnet, which excits the magnetic 
fields in the opposite direction. Therefore, the excited 
magnetic fields, B', are always anti-parallel in the space 
between neighboring silver layers and cancel each other. 
This explains the observation of a weak resonance with 
nearly zero n. For the symmetric mode shown in Fig. 
|6jb), the first and the fourth silver layers have current 
density along opposite directions and are almost uniform 
for all the metallic thickness of 30 nm silver layers. In 
the second and third silver layers the current density is 
no longer uniform in all thicknesses of the silver layers. 
Instead, the current flows along opposite directions on 
the two surfaces of each layer. Due to the anti-parallel 
current on the surfaces of the second and third silver lay- 
ers, the induced magnetic field, B', in the space between 
neighboring silver layers, is always parallel to each other. 
As a consequence, the 7 layers structure can be viewed as 
three cascade double-fishnet structures with the induced 
magnetic fields, B', along the same direction. Therefore, 
the symmetric mode results in a strong resonance with 
large negative n. Our detailed numerical work, shown in 
Fig. 4 and 6, explains very well why we obtain very low 
n ~0 at A=2230 nm and why we obtain high negative 
n = -2.5 at A=1859 nm. However, it is not clear that 
this current density distribution is responsible for high 
FOM shown in Fig. [5] 

The reason that the single unit cell (metal-dielectric- 
metal) has low FOM or high losses is due to its resonance 
structure. One way to increase the figure of merit, which 
is the ratio of |Re(n)/Im(n.)|, is to move away from the 
resonance frequency, where Im(n) is large, and therefore 
FOM can increase dramatically. This can be accoplished 
in both the weakly and strongly coupled fishnets, by in- 
troducing periodicity effects. For the single unit cell fish- 
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FIG. 7: The real parts of refractive index (red), Re(n), and 
the figure of merit (blue), for the single layer fishnet struc- 
tures with spacer thickness s = 0.025a (solid curves) and 0.1a 
(dashed curves), respectively. The other geometric param- 
eters are given by a x — a y — a, w x — 2a/5, w y — a/3, 
t = a/300, and the dielectric constant of the spacer is e r — 5. 

net, we can increase the size of the spacing layer and one 
can see from Fig. [7] that with a thicker spacing layer, 
s = 0.1a, the Re(n) reaches the Brillouin zone and the 
Re(n) < region is extend into the area where Im(n) ~ 0, 
so the FOM reaches a large value of 25.. 
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FIG. 8: The real (solid curves) and imaginary (dashed curves) 
parts of refractive index, Re(n) and Im(n), for the 3 (blue) 
and 19 (red) layers fishnet structures. The black dash line 
shows the position where Re(n) = — 1. 

In Fig. [51 we present both the real and the imagi- 
nary parts of the refractive index, Re(n) and Im(n), for 
the 3 and 19 layers fishnet structures. For the 3 lay- 
ers (the single layer of double-fishnet), the Re(n) has a 
smooth resonance curve (blue solid). The bandwidth of 



the Re(n) < region is relative narrow and close to the 
peak of Im(n) (blue dashed), so the figure of merit is very 
small (as shown in Fig. [5]). For the 19-layer fishnets, the 
Re(n) curve (red solid) does not have the resonance be- 
havior expected for a single functional layer, but it's very 
broad and has structure which is due to periodicity ef- 
fects [9]. Notice that for the 19 layer structure Re(«)=-1 
at A = 1688 nm and the Im(n) is 0.14, so the FOM is 
of the order of 10. However, for the 3 layer structure, 
Re(n)=-1 at A = 2075 nm and 2185 nm, and the lm(n) 
is 0.44 and 1.43 respectively, so the FOM is of the order 
of 1. Therefore, due to the distortion of Re(n) caused by 
the periodicity effects, the FOM of the fishnet structure 
increase dramatically as the number of layers increases. 



IV. CONCLUSIONS 

We have made a systematic study of the weakly and 
strongly coupled fishnets to understand the origin of neg- 
ative n, as well as the origin of losses and the large value 
of the FOM for the strongly coupled fishnets. We stud- 
ied the size dependence of the retrieved parameters (e, /i, 
and n) of the weakly and strongly coupled fishnet struc- 
tures. For both cases we found the retrieved parameters 
have a strong resonance behavior as the size of the unit 
cell decreases. We have also studied the convergence of 
the retrieved parameters, as the number of unit cells (lay- 
ers) increase. For the weakly coupled fishnet structures, 
we found the convergence results are relatively close to 
the single unit cell. Also, the converged FOM for the 
weakly coupled fishnet is the same order of magnitude as 
the single fishnet. For the strongly coupled fishnet struc- 
tures, we demonstrated that hybridization happens and 
we have two resonance modes. The antisymmetric reso- 
nance mode gives a strong negative n. As more unit cells 
or layers are added, the convergence of the retrieval pa- 
rameters are completely different than the single fishnet 
results and the FOM is much larger than the single fish- 
net. We have demonstrated that the large FOM for the 
strongly coupled fishnet is due to the periodicity effects. 
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APPENDIX A: EFFECTIVE PARAMETER 
RETRIEVAL FOR STRONGLY COUPLED 
FISHNETS 
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FIG. 9: Branches of the refractive index, Re(n) with m = 1 
(cyan), (blue), —1 (red), —2 (green), —3 (magenta) and —4 
(black). The cross, circle and diamond symbols in (a) rep- 
resent Re(n) for 7, 11 and 19 layer strongly coupled fishnet 
structure, respectively. The cross and circular symbols in (b) 
represent the 19 and 27 layer strongly coupled fishnet struc- 
ture, respectively. The shadow region shows where different 
branches overlap for 7, 11, 19, and 27 layers fishnet structures. 
The grey dotted lines show the branch boundaries which are 
given by m-K jkL. 

The effective retrieved parameter (e, /i, n and z) of 
single layer and many layers of metamaterial can be ob- 
tained from the transmission T and reflection coefficient 
R. There is a need for T and R to be inverted. As was 
discussed in detail in the literature [10-13], one can invert 
T and R 



z(uj) 
n(u>) 



= ±4 



'(l + i?) 2 -T 2 
(l-i?) 2 +T 2 



1 (l-R 2 +T 2 

±- — arccos — 

kL V 2T 



2tt 
'kL 



(Al) 
(A2) 



where L is the width of the homogeneous slab, and m = 
±1, ±2, ... . Note that both functions, z(ui) and n(uj), 
have multiple branches. The correct branch for z(lo) is 
chosen by imposing the physical requirement Re(z) > 
which is due to causality. The problem with the different 
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branches of Re(n) can be solved by considering different 
lengths for L, and one has to choose the branches that 
overlap. Especially if one has many layers, then many 
branches exist and one has to be very careful to select 
the correct ones. For the strongly coupled layers that the 
results were presented in Fig. 4, we would like to discuss 
how these branches were selected. The unit cell size is 
called do and it consists from metal-dielectric-metal and 
its width is do =160nm. 



1 




X(nm) 



FIG. 10: Transmission spectra for 7, 11, 19, and 27 layers 
strongly coupled fishnet system. 



In Fig. 9(a), we plot the branches and the retrieval re- 
sults for 7 layers (width= 2rfo), 11 layers (width= 3do) 
and 19 layers (width= 5do). Notice that the solutions for 
Re(n) overlap between 1200nm all the way to 2200nm 
and give negative values of Re(n). For A > 2200 nm, 
Re(n) ss and converges and the Im(n) w 3 in this 
region. So for A > 2200 nm the strongly coupled op- 
tical materials behave as a metal. In Fig. 9(b), we 
plot the branches and the retrieved result for 19 layers 
(width= 5do) and 27 layers (width= 7do) and one can 
see clearly that the convergence is much better for these 
larger systems. So we have solutions consisting of two dis- 
continued region for the Re(n), Re(rt) w for A > 2200 
nm and negative for 1500 nm< A < 2200 nm. 

In Fig. 10, we present the results for transmission, T, 
versus wavelength. Notice that for A > 2200 nm, T w 0, 
which is a metallic behavior and this is the reason that 
Re(n) k, and Im(n) « 3 for A > 2200 nm. 



